Decay of Affleck-Dine Condensates with 
Application to Q-balls 



Andrew Pawl a 

^Michigan Center for Theoretical Physics, Randall Laboratory, University of 
Michigan, 500 E. University Ave., Ann Arbor, MI 48109, USA 



Abstract 

Analytical and numerical estimates show that a charged Affleck-Dine condensate 
will fracture into Q-balls only when the Hubble time is significantly larger than 
the inverse soft-breaking mass of the field in question. This would generally imply 
that the decay of the field into light fermions will compete with Q-ball formation. 
We will show that for typical flat directions the large field value will significantly 
suppress decays of the condensate to fermions even if no baryon charge asymmetry 
exists. We will consider the details of the decay process for a condensate that does 
carry charge, and show that it is qualitatively different from that of an uncharged 
condensate. Finally, we will consider the possibility of resonant production of heavy 
bosons. We will show that this can have a strong effect on the condensate. Contrary 
to intuition, however, our results indicate that boson production would actually 
assist Q-ball formation in condensates with significant charge. 



1 Introduction 



A Q-ball is a non-perturbative solution of the equation of motion for a scalar 
field which is charged under a continuous U(l) symmetry [1]. The MSSM re- 
quires several such scalars, and so it is expected that Q-balls could be formed 
in a supersymmetric universe [2]. In fact, Q-balls seem inevitable in the con- 
text of Affleck-Dine [3] (AD) baryogenesis. Here, a flat direction composed 
of several squark fields gains a large expectation value and is set into coher- 
ent rotation by the action of a phase-dependent term in the potential. Such a 
charged scalar condensate has been shown to fracture into Q-balls in analytical 
treatments [4,5,6] and in numerical simulations [7]. 

Both numerical and analytical estimates agree that in gravity-mediated SUSY 
breaking models, Q-balls will only form when the age of the universe has 
reached a value of order ~ lC^m^ 1 where is the soft-breaking mass of 
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the AD condensate field (assumed to be of order 1 TeV). Naively speaking, 
however, we would estimate the decay width of the scalar field into light 
fermions to be of order: 



where we have substituted the gauge coupling g for the usual Yukawa coupling 
because the squark fields making up are coupled to gluino/chargino plus 
quark through gauge interactions. 

If we do not invent a suppression for this decay by making the fermions heavy 
then we see that the time for our condensate to decay can be as short as order 
lOOm^ 1 . Thus, we expect a competition between decay of the condensate into 
light fermions and fracturing of the condensate into Q-balls. As we will explore 
in this paper, however, there is not only a suppression to the decays, but 
there is also an important difference between the decay of a partially charge- 
asymmetric condensate and a neutral one. Taken together, these details will 
ensure that Q-ball formation is uninterrupted (and perhaps even aided) by 
decay of the condensate. 



2 Decay into Fermions 

2. 1 Necessity of Non-perturbative Approach 

The physics leading to the suppression of decay into fermions is familiar. Any 
fermions coupling directly to the field will gain a mass of order g\<f>\ where g 
is the relevant coupling and |0| is the magnitude of the complex scalar field. 
In modern formulations of the AD scenario [8], when the Hubble constant H 
reaches the AD field will have a magnitude of order: 



where M is a large mass scale (order Mqut or M p i) and n describes the flatness 
of the flat direction (two standard choices are n = 4 or n = 6). We will take 
n = 6 as our canonical value because it minimizes thermal concerns [9]. For 
typical numbers, then, we expect |0| > lO^m^ (or 10 9 m^, for n = 4) so that 
all fermions have effective masses much larger than the mass of the field. 

At first glance, this solves our problem completely since the field is stable. 
However, it is important that when H < the field will begin to execute 
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harmonic oscillations about = [8]. This will result in a sinusoidally vary- 
ing mass for the fields coupled to the field. Such a situation has already 
been studied for real (non-complex) scalar condensates in the context of post- 
inflation reheating, and has been shown to lead to decay of the condensate 
[10,11]. We wish to expand this analysis to the case of an oscillating complex 
field. 

Explicitly, we anticipate that after a few oscillations the field will see the 
effective potential: 

f/(|0|)^|0| 2 (3) 



where higher order terms can be neglected due to the small size of |0| [8]. 
We can now make the analogy between |0| and the radial position r of a 
particle in an r 2 potential. We know from basic classical mechanics that the 
angular momentum in such a system will be conserved (in this case, "angular 
momentum" is equivalent to baryon number), and further that the particle will 
follow closed orbits [13]. Numerical integration shows that this approximation 
is very nearly exact even in the presence of small corrections due to non- 
renormalizable terms and log running of the mass parameter. 

All of this amounts to the fact that we lose no generality in parameterizing 
the final solution for the field in the form of an ellipse centered on the origin: 

= asin(m^t) — ibcosim^t) (4) 



where we have assumed that a> b. 

Using the Noether current expression for the global U(l) yields the expression 
for net baryon number density of the condensate: 



n B 



i(3 (0*0 - 00*) = 2p m(f) ab (5) 



where (3 is the baryon charge per particle (usually 1/3). Then, by assuming 
that each scalar particle associated with the field has an energy of approxi- 
mately m^, it is simple to show that the ratio of net baryons to total scalars 
is: 

hb_ 2(3ba 

~ a 2 + b 2 ' U 

From this expression, we can see that the limit 6 = corresponds to a com- 
pletely uncharged scalar field. The limit b = a indicates a total charge asymme- 
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try (the condensate is made up entirely of baryons or entirely of antibaryons). 
Intermediate values of b indicate a partially charge-asymmetric condensate. 



Using the parameterization (4), we also find that: 

|0| = v /62 + ( a 2_ 6 2 )sin 2 (m ^ ) _ (7) 



Thus, we do anticipate an oscillating mass for any fields coupled to the AD 
scalar. 



2.2 Method of Imaginary Time 



The method of imaginary time is convenient for the calculation of non-perturbative 
production of fermions. We will use the results of [10] essentially verbatim. For 
an introduction to the imaginary time formalism, see the review [12] and the 
references therein. 

The method relies on finding the branch points in the complex-time plane of 
the fermion Hamiltonian. For our fermions, we expect: 

U = ( P 2 + 9 2 \<P\ 2 ) 1/2 = (p 2 + 9 2 b 2 + g\a 2 - b 2 ) sin 2 (m,t)) 1/2 . (8) 



Fortunately, exactly this form was treated in [10]. 



Two limits are analytically approachable using their results. First, note that 
by Equation (2) we always expect m$ <C a. Now, if we also assume b <^ a 
(equivalent to m <C mi in the conventions of [10]) we can take an analytical 
limit: 



e n / 2 m 3 J 2 
87r 2 {gay/ 2 



(9) 



(note that e here is the base of the natural log). In this approximation we have 
basically followed [12] exactly, except that we have added the assumption that 
log(4a/fo) is of the order n. Note that here we have assumed b is nonzero, and 
larger than the momentum p. If b were to approach zero, the results of [10,12] 
apply exactly (there is no divergence). Putting in typical numbers will tell us 
that T is of order 10~ 6 m^ at the largest (lO~ 5 m0 for n = 4). This suppression 
is enough to keep our condensate intact until Q-balls can form, even for the 
less-favorable n = 4 case. 
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Fig. 1. Exponential damping constant in the (j) — > tptp decay width for various values 
of b/a. ~ e _2< 3. We have assumed ga = 10 10 m^. 



Next let us take the limit b ~ a. Here again we can analytically approximate 
the decay width (this is equivalent to m mi in [10]): 



m 



3/2 
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which, for our typical numbers yields a decay width that has been exponen- 
tially suppressed to the extent that it is effectively zero. 

In each of these limits for the ratio of b to a, order of magnitude estimates 
show that the decay rate is highly suppressed. It is logical that the rate does 
not peak between, but rather goes smoothly from one to the other. To be safe, 
however, we can numerically evaluate the complete elliptical functions outside 
the range of validity of our approximations to estimate the exponential part 
of the suppression for arbitrary values of b/a. The curve shown in Figure 1 
computes only the exponential suppression of the decay rate, but it is enough 
to show that this assumption is correct. Thus, the large field values suppress 
the decay into fermions sufficiently to allow Q-ball formation to proceed. 



It is important to notice in Equation (10) that at b = a the decay rate goes to 
zero. This is no accident. A completely charge asymmetric scalar condensate 
is stable against decay into fermions - even non-perturbative decay [14]. This 
is an important component of Q-ball stability. 
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3 Parametric Resonance 



Finally, we must consider the decay of our field into bosons. This decay will 
be suppressed by the same physics that affected the decay into fermions. It is 
well known, however, that in the case of decay to bosons enhancement from 
particle statistics can offset this suppression. This phenomenon was mentioned 
in [10,11]. The first modern treatment, however, was [15]. We will follow the 
expanded version presented in [16]. 

We expect the field to have several couplings of the form g 2 0| 2 |x| 2 where g 
is the strong gauge coupling and x is a scalar (a linear combination of squarks 
that is orthogonal to the combination that makes up the direction). 

This contribution to the Lagrangian will result in an oscillating mass term for 
the x fidd when the <fi field is undergoing elliptical oscillation. The scalar x 
will have the following equation of motion (after transforming to momentum 
space): 



k 2 

— + g 2 b 2 + g 2 (a 2 - b 2 ) sin 2 (m^t) 



X = (11) 



where R is the scale factor of the universe (H = R/R). We can simplify this 
by using the substitution X = R?I 2 X- This gives: 



X + 



^ + g 2 b 2 + g\a 2 - b 2 ) sin \m,t) - \ H 2 - || 



X = 0. 



(12) 



Noting that the typical scale for H is <C a, b at the time 
oscillating, we assume we can rewrite this: 



is coherently 



X + 



R 2 



g 2 b 2 + g 2 (a 2 - b 2 ) sin 2 (m (/) t) 



X = 0. 



(13) 



This is almost exactly the case of stochastic resonance in the expanding uni- 
verse treated in [16]. The only issue to be careful of is the size of b. We will 
discuss that shortly. In the meantime, we will adapt section IX of [16] to our 
purposes. 

The fundamental procedure of the parametric resonance approach is the rewrit- 
ing of our equation of motion in the form of the Mathieu equation: 

x" + (A-2qcos(2z))x = (14) 
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where prime denotes differentiation with respect to z (please note that q in 
this context is an unfortunate choice for us - it has nothing to do with quarks). 
The appropriate substitutions here are: 



A = 




(15) 



Q = 



9 2 (a 2 -b 2 )_ 



(16) 



and 



z = mt. 



(17) 



The work of [16] shows that it is q which determines the efficiency of the con- 
densate decay. For q > 10 4 , we expect the condensate to retain a considerable 
q value even after completion of the "first stage" of preheating (with low den- 
sity of decay products). If q remains greater than order one after this phase is 
completed we expect the decay to continue and the condensate will give up a 
significant fraction of its energy to decay products. 

We are interested, therefore, in estimating q for a typical Q-ball scenario. At 
the beginning of oscillations we expect a 2 — b 2 to be of order lO^m^ which 
gives a q factor of: 



It is pointed out in [16] that the q = 10 4 cutoff value for strong preheating 
is weakly model dependent. Given that we expect g 2 ~ 0.1 for strong inter- 
actions, however, it seems safe to assume that we are well into the strong 
preheating regime (this is also true for the n = 4 flat directions, where we 
expect q ~ 10 14 g 2 ). 

Now let us consider, as we did for fermion production, what it would mean 
for this decay to go forward. Once again, we see that the decay efficiency is 
proportional to the difference a 2 — b 2 so that it vanishes as we approach a com- 
pletely asymmetric condensate. In fact, in the case of resonant production of 
bosons, we can have a second suppression since b is functioning as an effective 
bare mass in the formulas above. Thus, as mentioned in [16] we might expect 
preheating to become inefficient if 2b 2 > a 2 — b 2 . It appears, then, that our 
decays take away the neutral condensate but could leave a charged remnant. 

The physics behind this result is straightforward. The four point coupling 
fl ,2 |0| 2 |x| 2 that we have considered here should only mediate annihilation, not 



1n 20 2 

q ~ 11) g . 



(18) 
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true decays. This would certainly respect any baryon asymmetry present in 
the condensate. 

This leaves us with an important possibility. If decay into bosons is strong 
enough, it could be that generic AD condensates will damp much of their 
ellipticity. This process could actually aid Q-ball formation. In fact, even the 
time scales estimated in [16] for the decay process are right for Q-ball for- 
mation. They expect the first stage of resonant production to end at a time 
scale of order lOOm^ 1 , just before Q-balls would form according to numerical 
simulations [7]. 

There are two important issues to address, however, before we can say with 
confidence that annihilations help produce Q-balls. First, it is important to 
note that we have not included a study of rescattering of the produced x par- 
ticles. Such back-reaction could have a negative impact on Q-ball formation, 
and should be examined. Second, annihilation would almost certainly have a 
negative impact on Q-ball production in a very weakly charged condensate, 
since it could lead to significant decay of the condensate before Q-ball forma- 
tion would occur. 



4 Conclusion 



In summary, decay of a partly charge-asymmetric Affleck-Dine condensate 
into fermions will be strongly suppressed. This suppression will give more 
than enough time for the condensate to fracture into Q-balls. Annihilation 
of the neutral part of the condensate, however, can be enhanced by non- 
perturbative effects completely analogous to preheating in inflation. These ef- 
fects must be studied further, as they could have important consequences for 
Q-ball formation. In particular, it seems that they should make Q-ball forma- 
tion in strongly charged condensates more likely, while suppressing formation 
in weakly charged condensates. 
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